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Abstract 

The self-repelling Brownian polymer model (SRBP) initiated by Durrett and 
Rogers in [7] is the continuous space-time counterpart of the myopic (or 'true') 
self-avoiding walk model (MSAW) introduced in the physics literature by Amit, 
Parisi and Peliti in [1]. In both random motion in space is pushed towards 

domains less visited in the past by a kind of negative gradient of the occupation 
time measure. 

We investigate the asymptotic behaviour of SRBP in the non-recurrent dimen- 
sions. First, extending Id results from [21], we identify a natural stationary (in 
time) and ergodic distribution of the environment (essentially, smeared-out occupa- 
tion time measure of the process), as seen from the moving particle. As main result 
we prove that in three and more dimensions, in this stationary (and ergodic) regime, 
the displacement of the moving particle scales diffusively and its finite dimensional 
distributions converge to those of a Wiener process. This result settles part of the 
conjectures (based on non-rigorous renormalization group arguments) in [1]. 

The main tool is the non-reversible version of the Kipnis- Varadhan-type CLT 
for additive functionals of ergodic Markov processes and the graded sector condition 
of Sethuraman, Varadhan and Yau, [19]. 

MSC2010: 60K37, 60K40, 60F05, 60J55 
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1 Introduction and background 

The asymptotic scaling behaviour of self-repelling random motions with long memory 
has been a mathematical challenge since the early eighties. The two basic models con- 
sidered in the physical and probabilistic literature are the so-called myopic (or 'true') 
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self-avoiding random walk (MSAW), which appeared first in the physics literature in [1], 
and the so-called self-repelling Brownian polymer (SRBP) model, which was initiated in 
the probabilistic literature in [14], [7]. The two models (or, better said, families of mod- 
els), although having their origins in different cultures and having different motivations, 
are phenomenologically very similar. 

The simplest formulation of the MSAW model is as follows: Let X{n) be a nearest 
neighbour random walk on Z*^ and 

e{n, y) £(0, y) + | {0 < m < n : X{m) = y} \ (1) 

its occupation time measure, taken with some (possibly signed) initial values £(0, y) e Z. 
The walk is governed by the following law: 

P(X(n + l) = y I past, X(n) =x) = (2) 

r{({n..v)-({n.y)) 



1{ |a;-y| = l} 



^-x 1=1 ^) - ^)) 



where r : Z — > (0, oo) is a non-decreasing (and non-constant) weight function. (In the 
original [1] and the subsequent physics papers, the specific choice r{u) — exp{/3u}, /3 > 0, 
was made.) This means phenomenologically that the random walk X{n) is pushed by 
the negative gradient of its occupation time measure. 

The SRBP model is defined as follows. Let V : M*^ — > R+ be an approximate identity, 
that is a smooth (C°°), spherically symmetric function with sufficiently fast decay at 
infinity, and 

F:R'^^R'', F{x):^-grsAV{x). (3) 

For reasons which will be clarified later, we also impose the condition of positive definite- 
ness of V: 

V{p) := (27r)-''/' / e'P-''V{x)dx > 0. (4) 

A particular choice could be V{x) := exp{— | x f /2}. 

Let t ^ B{t) e M*^ be standard d-dimensional Brownian motion and define the 
stochastic process t ^ X{t) e R*^ as the solution of the SDE 



(5) 



X{t) = B{t) + [ [ F{X{s) - X{u)) du ds, 
Jo Jo 

or ^ 

dX{t) = dB{t) + ( [ F{X{t) - X{u)) du) dt. (6) 
Jo 

Remark: Other types of self-interaction functions F give rise to various different asymp- 
totics. For the few rigorous results (mostly in Id), see [14], [7], [3], [4] and in particular 
[13] which also contains a survey of the results. Recent Id results appear in [21]. 
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Now, introducing the occupation time measure 

£{t, A) := ^) + I {0 < s < t : X{s) G A} \ (7) 

where A C M*^ is any measurable domain, and i{0, A) is some signed initialization, we 
can rewrite the SDE (6) as follows: 

dX{t) = dB{t) - grad {V * £{1, ■)) {X{t)) dt (8) 

where * stands for convolution in R''. We assume that i{0,A) is a signed Borel measure 
on with slow increase: for any s > 

lim N-^'^+'^ \£\{0, [-N, Nf) = 

Af^oo 

The form (8), compared with (2), shows explicitly the phenomenological similarity of the 
two models. 

Non-rigorous (but nevertheless convincing) scaling and renormalization group argu- 
ments originally formulated for the MSAW, but equally well applicable to the SRBP 
suggest the following dimension-dependent asymptotic scaling behaviour (see e.g. [1], 
[15], [17]): 

- in d — 1: X{t) ~ t"^^^, with intricate, non-Gausssian scaling limit; 

- in (i = 2: X(t) ~ t^/^(logt)'', with some controversy about the value of the exponent ( 
in the logarithmic correction, and Gaussian (that is Wiener) scaling limit expected; 

- in 0? > 3: X(t) ~ i^^^, with Gaussian (i.e. Wiener) scaling limit expected. 

Ind — 1, for some particular cases of the MSAW model (MSAW with edge, rather than 
site repulsion and MSAW with continuous time and site repulsion), the limit theorem 
for X(n)/n^/'^ was established in [23], respectively, [24], with the truly intricate limiting 
distribution identified. The scaling limit of the process t i— > X{nt)/n^/'^ was constructed 
and analyzed in [25]. The proofs in [23] and [24] have some built-in combinatorial elements 
which make it difficult (if possible at all) to extend these proofs robustly to a full class 
of Id models of random motions pushed by the negative gradient of their occupation 
time measure. However, more recently, a robust proof was given for the super-diffusive 
behaviour of the Id models: in [21], inter alia, it is proved that for the Id SRBP models 
\\m^^^t-^'^^{X{tf) > 0, nmt^oo^"^/^E(X(t)2) < oo. These are robust super-diffusive 
bounds (not depending on microscopic details) but still far from the expected t'^^'^ scaling. 

In d — 2, very little is proved rigorously. For a modified version of MSAW where self- 
repulsion acts only in one spatial (say, the horizontal) direction, the marginally super- 
diffusive lower bound lim^_,,, t-H\oKt)~^/^B(X(tY) > holds, cf. [26]. 

In the present paper, we address the d > 3 case of the SRBP model. We identify a 
stationary and ergodic distribution of the environment as seen from the position of the 
moving point and in this particular stationary regime, we prove diffusive limit (that is 
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non-degenerate CLT with normal scaling) for the displacement. Our general approach is 
that of martingale approximation for additive functionals of ergodic Markov processes, 
initiated for reversible processes in the classic Kipnis - Varadhan paper [10] and extended 
to non-reversible cases in [22], [27], [19]. We shall refer to this approach as the Kipnis - 
Varadhan theory. In particular, validity of the efficient martingale approximation will 
rely on checking the graded sector condition of [19]. 

Similar results for the MS AW model on the lattice Z"^, (i > 3, will be presented in [8] 
Next, we describe our results in plain words. For precise formulations, see subsection 

2.2. 

As a first step, we note that the environment profile appearing on the right-hand 
side of (6), (8), as seen in a moving coordinate frame tied to the current position of the 
process, 1 1— > r){t, •): 

rj{t,x) :^r]{0,X{t)+x)+ [ V{X{t) + x - X{u)) du (10) 

Jo 

= 77(0, X{t) +x) + {V* i{t, ■)) {X{t) + x) 

is a Markov process in a properly chosen function space Q, to be specified later. As a first 
step, we identify a natural time-stationary and ergodic distribution of this process (10). 
Rather surprisingly, this is the Gaussian (scalar) field x h- > uj{x) e M with expectation 
and covariances 

E(^(x)) = 0, C{x - y) := ^{uj{x)uj{y)) =g*V{x-y) (11) 

where 

^:R''^M, g{x) ■.^\x\^~'^ (12) 

is the Green function of the Laplacian in W^. Note that throughout this paper d > 3. 
This is the massless free Gaussian field whose ultraviolet singularity is smeared out by 
convolution with the smooth and rapidly decaying approximate identity V . The Fourier 
transform of the covariance is 

C{p) = \p\-'^V{p). (13) 

See Theorem 1. All further results will be meant for the process being in this stationary 
regime. From this result, by ergodicity, the law of large numbers for the process X{t) 
drops out, see Corollary 1. 

The main result of the paper refers to the diffusive limit of the process t i— > X{t). 
From (5) and (10), it arises that the displacement is written as 

X{t) = B{t)+ [\{v{s))ds (14) 

where 99 : — M"' is a function of the state of the stationary and ergodic Markov process 
t ^ r]{t): 

cp{u;) = -gradu;(0). (15) 
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So, the natural approach to the diffusive limit of X{t) is the Kipnis - Varadhan theory. 
We will prove validity of an efficient martingale approximation by checking Sethuraman- 
Varadhan- Yau's graded sector condition, cf. [19]. It is easy to see that due to the 
spherical symmetry of the problem, we get 

E{Xkit)Xi{t)) = 5k,id-'E{\X{t) n. (16) 

We prove that for d> 3, the limiting variance 

(7^ := d-^ lim r^E( I X{t) 1^ ) e (0, oo) (17) 

t— >oo 

exists and the finite dimensional marginals of the diffusively rescaled process 

X(Nt) , , 

Xnit) := (18) 
ay N 

converge to those of a standard d-dimensional Brownian motion. See Theorem 2. The 
main result shows similarity in spirit and techniques with those of [12], but the differences 
are also clear. 

The results are meant in probability with respect to the initial profile rj{0,x) sampled 
from the stationary (and ergodic) initial distribution hinted at above. Recent results 
by Cuny and Peligrad [5] raise the hope that the Kipnis - Varadhan theory could be 
enhanced to CLT for almost all initial conditions sampled according to the stationary 
distribution. 

The rest of the paper is structured as follows: in section 2, we give the formal def- 
initions, introduce notations, identify the stationary measure and formulate our main 
results precisely. In section 3, we give the functional analytic background: the Hilbert 
spaces and the (bounded and unbounded) linear operators involved and the infinitesi- 
mal generator will be presented. Ergodicity and LLN for the displacement drop out for 
free. The short section 4 is devoted to recalling the martingale approximation in Kipnis - 
Varadhan theory and the graded sector condition of [19]. Finally, in section 5, we check 
the abstract functional analytic conditions for our particular problem, and we conclude 
with the proof of the CLT for the displacement. 

2 Formal setup and results 
2.1 The stationary measure 

We start with the Ansatz that the stationary distribution of the process t i— > ri{t, ■) is 
translation invariant zero mean Gaussian scalar with some covariance 

E{rj{t,x)r]{t,y))=C{y-x), (19) 
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to be identified at the end of tlie following computations. 

In order to prove that this is indeed time-stationary, we have to show that for any 
test function x — > u{x) £ M, the moment generating functional 

(t>{t,u) :=E(exp{(ii,77(t))}) (20) 

is actually constant in time. In the present subsection, we use the notation 

{u,v) :— I v{x)u{x)dLX. (21) 

In the forthcoming computations of the present section all repeated subscripts are 
summed from 1 to d. Using (10), note that, by standard Ito calculus, 

d(«, r]{t)) = -{dm, r]{t)) dBi{t) + ^{dfiu, r]{t)) dt - {dm, v{t))div{t, 0) dt + {u, V) dt. 

(22) 

Hence 

E{deM{%vit))}\^t)^ (23) 
exp{(«,7y(i))} (^l{^iu,ri{t)) + l{dm:V{t)f - {diu,ri{t))diri{t,0) + {u,V)^ dt. 

Now, using the Ansatz that x ^ r){t, x) (with t fixed!) is a Gaussian field with covariance 
(19), by standard computations of Gaussian expectations, from (23), we obtain 

dE{eM{uMt))}) ^ ^^p^^^^ ^ ^ ^^/2} [^{dlu, C*u) + \ {dm, c * dm)+ (24) 

+ ^{dm, c*u)^- {dm, diC) - {dm, c * u) {u, dfi) + {u, v)) dt. 

On the right-hand side the first two terms cancel out by an integration by parts. The 
third and fifth terms cancel one by one due to the simple fact that for any test function 
u, 

{dm, C * -u) = 0. (25) 
Thus, the right-hand side of (24) is canceled out completely iff 

V = -die. (26) 

This is equivalent to (11), (13). Note that o? > 3 was assumed. 

2.2 Formal setup and results 
2.2.1 State space and Gaussian measure 

The proper state space of our basic processes will be the space of smooth scalar fields of 
slow increase at infinity: 

Q := {a; e C°°(W^ ^ R) : \\uj ||^^^ < oo} (27) 



where 11^*^11^^ are the seminorms 

-l/r 



I Wm.T 



sup {\^\x\y"' \d\Zl,n.A^)\ (28) 



defined for the multiindices m = (mi, . . . , m^), rrij > 0; and r > 1. The space fl endowed 
with these seminorms is a Frechet space. 

From Minlos's theorem (see [20]), it follows that there exists a unique Gaussian prob- 
ability measure 7r(da;) on the space of tempered distributions S'{W^ — > M) with charac- 
teristic functional 

E(exp{i(M,a;)}) = expi / / u{x)C{x - y)u{y) dxdyl , (29) 

and from smoothness of the covariance function C{x), it follows that the probability 
measure 7r(da;) is actually concentrated on the space ft C <S'(M'^ — > M) and (11) holds. 
The Gaussian field u!{x) is realized e.g. as a moving average of white noise: 

u;{x) = [ U{x-y)w{y)dy, (30) 

where U is the unique positive definite function for which U * U = V and w is d- 
dimensional white noise. 

The group of spatial translations 

3 z^T^-.n^n, {t^uj){x) := uj{x + z) (31) 

acts naturally on Q and preserves the probability measure 7r(da;). Actually, the dynamical 
system (Q, 7r(da;), : z e W^) is ergodic. 

2.2.2 Processes 

First, we consider the process 1 1— > {X{t), (^{t, •)) e R*^ x Q defined as follows: 

X{t)^X{0) + B{t)- [ gradC(s,X(s))ds, (32) 

Jo 

at, x) = C(0, x) + f V{x - X{s)) ds (33) 
Jo 

where t ^ B{t) is a standard rf-dimensional Brownian motion, and ^(0) £ W^, ^(0, •) G O 
are the initial data for the process t ^ (X(t), ^(t, ■)). 

Written as a single SDE for the process t ^ X{t), we get from (32), (33): 



X{t) = X(0) + Bit) + |c(0, X{s)) + F{X{s) - X{u)) dwj ds. 



(34) 



7 



The SDE (34) differs from the original SDE (5) only by the presence of the initial profile 
((0,x), which is a natural modification. 
From (32) and (33), it follows that 

X{to + t)=X{to) + {B{t + to)-B{to))- gradC(s,X(s))ds, (35) 

rta+t 

C{to + t,x)^C{to:X)+ V{x-X{s))ds. (36) 

From this form, it is apparent that the process t {X{t), ({t, ■)) G M*^ x f2 is Markovian. 
The environment profile as seen from the moving point X{t) is 

x^rj{t,x) :^C{t:X{t)+x). (37) 

From (35), (36), we readily obtain that t ri{t) := ri{t, ■) is itself a Markov process on 
the state space Q. 

We define the function (/? : Q ^ R'^ by (15), and from (32), (33), (37), we readily get 
(14). 



2.2.3 Results 

Theorem 1. The Gaussian probability measure -niduj) on Vl, with mean and covariances 
(11) is time-invariant and ergodic for the fl-valued Markov process 1 1— >• r]{t). 

Corollary 1. For n -almost all initial profiles ^(0, •), 

, X{t) , , 

lim = a.s. (38) 

t— »oo t 



Remarks: It is clear that, in dimensions ci > 3, other stationary distributions of the 
process t i— > r]{t) exist. In particular, due to transience of the process t h- > X{t), the 
stationary measure (presumably) reached from starting with "empty" initial conditions 
77(0, x) = certainly differs from our d7r. Our methods and results are valid for the 
particular stationary distribution d7r. 

The main result of the present paper is the following theorem: 

Theorem 2. In dimensions d> 3, the following hold: 
(i) The limiting variance 

(7^ := d-^ lim r^E( I X{t) f ) (39) 

exists and 

1 < CT^ < 1 + p2 (40) 
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where 

p'^-.^d-^ |p|-2y(p)dp< oo. (41) 
(a) The finite dimensional marginal distributions of the diffusively res caled process 

X.it) :^ ^ (42) 



converge to those of a standard d-dimensional Brownian motion. The convergence is 
meant in probability with respect to the starting state ri{0) sampled according to dvr. 

Theorem 2 will be proved by use of the martingale approximation of the Kipnis- 
Varadhan theory and the so-called graded sector condition of [19]. 

3 Spaces and operators 

The natural formalism for the proofs of our theorems is that of Fock spaces and Gaussian 
Hilbert spaces, and linear operators over them. For basics of Gaussian Hilbert spaces and 
Wick products, see [9], [20]. Our main Hilbert space is H := C'^{Q, n). This is a Gaussian 
Hilbert space, and has very natural unitary equivalent representations as Fock spaces. We 
follow the usual notation of Euclidean quantum field theory, see e.g. [20]. In subsection 
3.1, we give formal definition of the three unitary equivalent representations of the Hilbert 
space £^(fi,7r). In subsection 3.2, we define the linear operators which are relevant for 
our purposes and we present their action on the three imitary equivalent formulations. 
In subsection 3.3, the infinitesimal generator of the semigroup of the stationary Markov 
process t ^ vi'^r) ^ ^5 acting on £^(f2,7r), and its adjoint is computed and the first 
consequences (ergodicity, LLN) are settled. 

3.1 Spaces 

Throughout this paper, we use the convention of unitary Fourier transform 

u{p):^{2n)-'^/^ [ e'P-'^uix) dx. (43) 

and the shorthand notation 

d 

X = (Xi, ...,Xn)e {R'^y, Xm = (Xml, ■ ■ ■ , Xmd) £ ^'^ , d^l := ^ , (44) 

dXml 

p = (Pl,...,p„) G (M'^)", Pm = (Pml,---,Pmd) e K'^, (45) 

m = 1, . . . , n, / = 1, . . . , 

We denote by respectively, 5„, the symmetric Schwartz spaces 

Sn : M'^'" ^ C : 'u(otx) = ^(x), w G Perm(n)}, (46) 
Sn -.^{u : R''" ^ C : w(rop) = u{^), w e Perm(n)}. (47) 
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In the preceding formulas Perm(n) denotes the group of permutations on the n indices. 
The spaces »S„, respectively, »S„ are endowed with the following scalar products 

{u,v):=[ [ ^^C(x - y)^;(y) dxdy, (48) 
{u,v) := / ^a(p)^;(p)dp (49) 

n n 

C(x - y) := n ^(^- - ^-)' ^(P) •= n ^^Prn). (50) 



where 



m=l m=l 



Let Kn and Kn be the closures of respectively, Sn with respect to the Euclidean 
norms defined by these inner products. The Fourier transform (43) realizes an isometric 
isomorphism between the Hilbert spaces /C„ and /C„. 

These Hilbert spaces are actually the symmetrized n-fold tensor products 

:= symm(K:r) , X:„ := symm(£r) . (51) 

Finally, the full Fock spaces are 



The Hilbert space of our true interest is 7^ = tt). This is itself a graded Gaussian 

Hilbert space 

n = ®~o^n (53) 

where the subspaces Tin are isometrically isomorphic with the subspaces /C„ of /C through 
the identification 

(pn-JCn^Tin, 4>„{u) := ^= [ u{:>i) :uj{xi) . . . uj{xn): dx. (54) 

V n\ J]S.dn 

Here and in the rest of this paper, we denote by :Xi . . the Wick product of the 

jointly Gaussian random variables (Xi, . . . ,Xn). In order to ease notation the mapping 
01 : /Ci ^ Hi will be simply denoted by 0. 
As the graded Hilbert spaces 



Ti. '■— ©^o"^"' ^ '~ ^nLo^m ^ '■— ®^o^" (55) 

are isometrically isomorphic in a natural way, we shall move freely between the various 
representations. 



10 



3.2 Operators 
3.2.1 General notation 

We use the standard notation of Fock spaces. First, we give a general framework of 
notation and identities formulated over the Gaussian Hilbert space Ti. Then, we turn to 
our relevant linear operators and we give their representations in all three Hilbert spaces 
7i, /C and K. 

The action of linear operators over Ti = (B'^=o'Hn will be typically given in terms of 
Wick monomials. It is understood that their action is extended by linearity and graph 
closure. 

Given a (bounded or unbounded) densely defined and closed linear operator A over 
the basic Hilbert space /Ci, its second quantized version, acting over the graded Gaussian 
Hilbert space H will be denoted by dr{A). This latter one acts over Wick monomials as 

follows, dr(^) : Hn n„, 

n 
m=l 

Given a vector u from the basic Hilbert space ICi, the creation and annihilation 
(raising and lowering) operators associated to it, acting over the Gaussian Hilbert space 
H, will be denoted by a*{u) : Hn — Hn+i, respectively, a{u) : Tin — Ti-n-i, acting on 
Wick monomials as 

a*{u) -.(pivi) . . .(p{vn): = ■.(p{u)(p{vi) . . . (p{vn): , (57) 

n 

a{u) :(f){vi) . . . (f){vn): = ^{u, Vm) :0(^^i) ■ ■ ■ (f){Vm-l)(f>{Vm+l) ■ ■ ■ (Pivn)- ■ (58) 

m=l 

For basics about creation, annihilation and second quantized operators, see e.g. [20] or 
[9]. 

We also define the unitary involution J on 7i: 

Jf{uj):=f{-uj), J\n„=hirnn„. (59) 
The well-known canonical commutation relations between the operators introduced 

are: 

[a{u), a{v)] = 0, [a*{u), a*{v)] = 0, [a{u), a*{v)] = {u, v)I, (60) 

[dT{A), a*{u)] = a*{Au), [dT{A), a{u)] = ~a{A*u), (61) 

[J, dr{A)] = 0, {J, a*{u)} = 0, {J, a{u)} = 0. (62) 

Two more operators will be needed: given an element u E ICi, multiplication by (j){u) 
will be denoted M{u), that is, formally, for / e >C^(r2,7r), 

{M(u)f)(u;):^cl^{u){u;)f(u;). (63) 



11 



Finally, for a fixed element e Jl, we introduce differentiation in the direction formally 

D^f{uj) := lim£-i(/(a; + e^) - f{uj)). (64) 



£^0 

Both operators are well-defined on Wick monomials, and are extended by linearity and 
graph closure. 

Given u G /Ci the identities (65) and (66) below hold: 

(1) The multiplication operator M{u) is actually 

M{u)^a*{u) + a{u). (65) 

(2) liC*ueVL then 

Dc*u = a{u). (66) 

Both identities are checked by direct computation on Wick monomials. The identity (66) 
is a particular case of the directional derivative of Malliavin calculus, see [9]. 

3.2.2 Specific linear operators 

The most relevant operators for our present purposes are 

d 

Vi:^dr{di), A:=J]Vf, := a(9,5o), :^ a*{diSo) (67) 

1=1 

where di = -£r^ and Sq is Dirac's delta concentrated on e W^. Note that Sq and all its 
partial derivatives are in the Hilbert space /Ci. 

We give now their action on the spaces 7i„, /C„ and /C„. The point is that we are 
interested primarily in their action on the space £^(f2,7r) = ©^g^n, but explicit com- 
putations in later sections are handy in the unitary equivalent representations over the 

space /C = ©^q/C„. The action of various operators over Tin will be given in terms of the 
Wick monomials :a;(xi) . . .a;(xn): and it is understood that the operators are extended 
by linearity and graph closure. 

• The operators V;, I — 1, ... ,d: 

n 

'^I'-Ti-n^ Tin, V/ :a;(xi) . . . uj{xn): = - ^ :a;(xi) . . . diuj{xm) ■ ■ ■ uj{xn): , (68) 

m=l 

Wi-.Kn^lCn, V;m(x) = ^-^(x), (69) 

m=l 

n 

Vi-.lCn^iCn, Viu{p) = ^ ( ^ Pml) u{p) ■ (70) 

m=l 

Note that these are actually unbounded, closed, skew self-adjoint operators. They 
are densely defined on Ttn, A^n, respectively, /C„. 
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• The operator A: 



A-.Tin^Hn, Aw{xi)...u;{xn):= (71) 

d n 

^ ■.Uj{xi)...diUj{Xm)---diUj{Xm')---Uj{Xn)-, 



1=1 m,m'=l 

d n 



d u 



I , dx„adxm'i 

1=1 m.m =L 



(x), 



A : /C„ ^ /C„, Am(p) 



n 
m=l 



m(p). 



(72) 
(73) 



The operator A is unbounded, densely defined, self-adjoint and positive. Note that 
A is not the second quantized Laplacian. 



• The operator 


1 A|" 


-1/2 ^ 


(-A) 


1 A 


-1/2 


: Hn — 




1 A 


-1/2 


■ — 




1 A 


-1/2 


■ ICn — 





no explicit formula, 
no explicit formula. 



A 1-1/2 ^/ 

A I ' m(p 



n 
m=l 



(74) 
(75) 

u{p). (76) 



1/2 

The operator | A | 'is unbounded, densely defined, self-adjoint and positive. 
The operators | A |~^^^ V;, I — 1, ... ,d: 



1 A 






no explicit formula. 


1 A 






no explicit formula. 


1 A 






\A\-'/'Vl^p) = '^^'^"'\u{p) 

1 2^m=l PiT^ 


These are 


hounded skew self-adjoint operators with operator norm 



(77) 
(78) 

(79) 



(80) 
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• The creation operators a* , I — 1, . . . ,d: 

a*i : Hn T-Cn+i, a*i -.ulxi) . . . uj{xn): =:diu;{0)u;{xi) . . . a;(a;„): , 



(h ■ JCn /Cn+1, aiU{xi, Xn+l) = 



(81) 
(82) 



^ n+l 

, ^ y] did{xm)u{xi 

y/n + 1 

* m=l 



, . . . , Xjn—i, Xjn-^i, . . . , Xn-\-l)j 



Gi : ICn ^ }Cn+l, a*iU{pi, . . . ,Pn+l) 



(83) 



^ n+l 

== X] iPmlu{Pl, ■ ■ ■,Pm-l,Pm+l, ■ ■ ■ , Pn+l) ■ 

m=l 



The creation operators a^*, restricted to the subspaces 7Y„, /C„, respectively, /C„ are 
bounded, with operator norm 



= ^/c(o)V^H^. 



(84) 



• The annihilation operators o/, / = 1, . . . ,d: 

ar.Hn^ Hn-i, ai :uj{xi) . . . Lu{xn): = (85) 

n 

^ diC{xm) -.ivixi) . . .u;{xm-i)uj{xm+i) ■ ■ ■^^(a^n): , 

m=l 

ai : fCn ^ ICn-i, aiu{xi,...,Xn-i) ^ Vn u{xi, . . . ,Xn-i,y)diC{y) dy, (86) 

jRd, 

ai : ]Cn ^ JCn^i, aiu{pi,...,Pn-i)^Vn u{pi, . . . ,Pn-i,q)iqid{q) dq. (87) 



The annihilation operators ai restricted to the subspaces Tin, K.n, respectively, /Cj, 
are bounded with operator norm 



ai 



V^CXO)^/^. 



(88) 



Furthermore, as the notation and ai suggests, these operators are adjoint of each 
other. 



Since all computations will be performed in the representation /C, we give a common 
core for all the unbounded operators defined above - and some others to appear in future 
sections: 

C := ©^=o^n, Cn := {m e : sup | m(p) | < oo}. (89) 
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1/2 

Note that the operator | A | ' is defined on the dense subspace C only for d > 3. 
Furthermore, in dimensions d > 3, the operators | A |~^^^ defined on the dense 

subspaces Cn, are essentially self-adjoint. This follows, e.g., from Propositions VIII. 1, 
VIII.2 of [18]. 

Notice also that V is the infinitesimal generator of the unitary group of spatial trans- 
lations while A is the infinitesimal generator of the Markovian semigroup of diffusion in 
random scenery 

exp{^V} = r„ r,/(^) := /(r,a;), (90) 

exp{tA} = g„ Qtf{uj) := / exp{-^V(2t)} ^^^^^^ ^^^^ 

3.3 The infinitesimal generator, stationarity, 
Yaglom reversibility, ergodicity 

We denote by Pt the semigroup of the process r]{t): 

P^:n^n, Ptfiu) := E{f{v{t)) I v{0) = cv). (92) 

Then [0, oo) 3 t Pt E BiTi) is a Markovian contraction semigroup on Ti. In order to 
identify its infinitesimal generator, note that the infinitesimal change in the state of the 
Markov process r]{t) is due to the following three terms: 

(1) infinitesimal spatial shift due to dB{t); 

(2) infinitesimal spatial shift due to — grad77(i, 0)di; 

(3) infinitesimal local change in r] due to increase of local time. 

Altogether 

r]{t + dt, x)^r]{t,x + dB{t) - grad r]{t, 0)dt) + V{x)dt. (93) 
Hence, given a sufficiently regular function on the state space / : Q — > R, we compute 

lim^ ^ ■ ^ = I -A - 2_^M(9i5o)Vi + Dy I /(u;). (94) 

Recall (11) and note that hence 

d 

V^-C*Y.^l5o (95) 



with 



-Y^dldoelC^. (96) 



1=1 



15 



Using (65), (61) and (66) (in this order), we readily obtain the following expression for 
the infinitesimal generator of the semigroup P^: 

1 

G':=-A + J](a;Vi + V^aO. (97) 
1=1 

This operator is well defined on Wick polynomials of the field ui{x) and is extended by 
linearity and graph closure. It is not difficult to see that it satisfies the criteria of the 
Hille - Yoshida theorem (see [18]) and thus it is indeed the infinitesimal generator of a 
Markovian semigroup. We omit these technical details. 
The adjoint generator is 

1 

G* -.^ -A-Y,{<^i + ^iai). (98) 

1=1 

Note that due to the inner coherence of the model the last two terms on the right hand 
side of (94) combine to give the tidy skew self-adjoint part of the infinitesimal generators 
in (97), (98). 

For later use, we introduce notation for the symmetric (self-adjoint) and anti-symmetric 
(skew self-adjoint) parts of the generator 

5:^_1(G + G'*) = -1a, (99) 

1 

-(G'-G'*) = ^(a;V, + Via,)=:^+ + ^_. (100) 
1=1 

It is a standard - though not completely trivial - exercise to check that the operators 
S and A, a priori defined on the dense subspace C are indeed essentially self-adjoint, 
respectively, essentially skew self-adjoint. 
Note that 

S :Hn^ 'Hn: Aj^ : Hn ^ "^n+l, ^- '■ "^n ^ "^n-l, ^=F — ~^±> (101) 

and 

S \no= 0, ^^0= 0, A_ \ho®h,= 0. (102) 

It is clear that 

G*! = 0, (103) 

and hence, it follows that tt is indeed stationary distribution of the process t ^ ri{t) and 
G* is itself the infinitesimal generator of the stochastic semigroup P^ of the time-reversed 
process. 

Actually, so-called Yaglom reversibility holds. From (62), it follows that 

G* = JGJ. (104) 
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This identity means that the stationary forward process (—00, 00) 3 t ^ rjlt) and the 
flipped backward process 

(-00, 00) 9 t ^ rj{t) := -vi-t) (105) 

obey the same law. This is a special kind of time-reversal symmetry called Yaglom 
reversibility, see [28], [29], [6]. 

Proving ergodicity is easy: The Dirichlet form of the process t 1— > r]{t) is 



V(f) -(/, Gf) = A/) = ^ II V J 11^ . (106) 

1=1 

So, 

{P(/)=0} « {VJ = 0, / = 1,...,/} ^ {/ = const. TT-a.s.}, (107) 

since z ^ acts ergodically on {Q, tt). 

This proves Theorem 1. Corollary 1 follows directly from (14) by the ergodic theorem. 



4 CLT for additive functionals of ergodic Markov pro- 
cesses, graded sector condition 

In the present short section we recall the non-reversible version of the Kipnis - Varadhan 
CLT for additive functionals of ergodic Markov processes and the graded sector condition 
of Sethuraman, Varadhan and Yau, [19]. 

Let (fi, J^, tt) be a probability space: the state space of a stationary and ergodic 
Markov process t 1— ri{t). We put ourselves in the Hilbert space H := C^{Q,n). Denote 
the infinitesimal generator of the semigroup of the process by G, which is a well-defined 
(possibly unbounded) closed linear operator on Ti. The adjoint generator G* is the 
infinitesimal generator of the semigroup of the reversed (also stationary and ergodic) 
process ri*{t) = rj{—t). It is assumed that G and G* have a common core of definition 
C <^ Tt. Let f & H, such that (/, 1) — f dir — 0. We ask about CLT/invariance 
principle for 

N-"^ fHs))ds (108) 
^0 

as N ^ 00. 

We denote the symmetric and anti-symmetric parts of the generators G, G*, by 

S ■= -hG + G*), A:=Ug-G*). (109) 

These operators are also extended from C by graph closure and it is assumed that they 
are well-defined self-adjoint, respectively, skew self-adjoint operators 

S*^S>0, A*^ -A. (110) 



17 



Note that —S is itself the infinitesimal generator of a Markov semigroup on tt), for 

which the probability measure tt is reversible (not just stationary). We assume that —S 
is itself ergodic: 

Kct{S) = {cl : ceC}. (Ill) 
We denote by Rx e B{7i) the resolvent of the semigroup s i— > e**^: 

poo 

Rx -.^ e-^'e'^ds ^ {XI - G)~\ A > 0, (112) 
Jo 

and given / e 7i as above, we will use the notation 

ux-Rxf. (113) 

The following theorem yields the efficient martingale approximation of the additive 
functional (108): 

Theorem (K V) . With the notation and assumptions as before, if the following two limits 
hold in Ti, 

lim A^/^-ua = 0, (114) 



A-»0 

then 



\imS^^\x =■■ V e n, (115) 
(7^=21im(iiA,/) e [0,oo), (116) 

A — >'0 

and there exists a zero mean, JC^ -martingale M(t), adapted to the filtration of the Markov 
process r](t) with stationary and ergodic increments and variance 

E(M(t)2) = aH (117) 

such that 

lim Ar-iE(( / f{r]{s))ds- M{N)Y) ^Q. (118) 



In particular, if a > 0, then the finite dimensional marginal distributions of the rescaled 
process t i— > a'^N'^^^ f{Tj{s)) ds converge to those of a standard Id Brownian motion. 

Remarks: 

(1) Conditions (114) and (115) of the theorem are jointly equivalent to the following 

lim (A + A')(^^A,^iA') = 0. (119) 

A,A'^0 

Indeed, straightforward computations yield: 

(A + >!){ux, ux') = II S^/\ux - ux') f + \\\ux f + A' || uy f . (120) 



18 



(2) The theorem is a generalization to non-reversible setup of the celebrated Kipnis- 
Varadhan theorem, [10]. To the best of our knowledge, the non-reversible formulation, 
proved with resolvent rather than spectral calculus, appears first - in discrete-time 
Markov chain, rather than continuous-time Markov process setup and with condition 
(119) - in [22] where it was applied, with bare hand computations, to obtain CLT for a 
particular random walk in random environment. Its proof follows the original proof of 
the Kipnis - Varadhan theorem with the difference that spectral calculus is to be replaced 
by resolvent calculus. 

(3) In continuous-time Markov process setup, it was formulated in [27] and applied to 
tagged particle motion in non-reversible zero mean exclusion processes. In this paper, 
the (strong) sector condition was formulated, which, together with an if_i-bound on the 
function / G 7Y, provide sufficient condition for (114) and (115) of Theorem KV to hold. 

(4) In [19], the so-called graded sector condition is formulated and Theorem KV is 
applied to tagged particle diffusion in general (non-zero mean) non-reversible exclusion 
processes, in d > 3. 

(5) For a more complete list of applications of Theorem KV together with the strong 
and graded sector conditions, see the surveys [16], [11]. 

Checking conditions (114) and (115) (or, equivalently, condition (119)) in particular 
applications is typically not easy. In the applications to RWRE in [22], the conditions 
were checked by some tricky bare hand computations. In [27], respectively, [19], the 
so-called sector condition, respectively, the graded sector condition were introduced and 
checked for the respective models. 

We recall from [19] the graded sector condition. Assume that the Hilbert space 

n = tt) is graded 

n = ©~ o^n, (121) 

and the infinitesimal generator is consistent with this grading in the sense of (101). 

Theorem (SVY). Assume that the Hilbert space and the infinitesimal generator G = 
—S + A are graded in the sense specified above and, in addition, there exist G [0, 1) and 
C < oo such that for any n E N and any g G Tin, h G Ti-n+i 

I {h,A+g) I < Cn^./{h;Sh:j./{^). (122) 

If f &TI with (/, 1) — is such that 

11^-1/2/ II :=lim(/,M,)<oo, (123) 

then (114) and (115) hold and, as consequence, the conclusions of Theorem KV are valid. 
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5 Proof of the CLT 



The proof of Theorem 2 consists of three parts. In paragraph 5.1.1, we prove diffusive 
lower bound on the variance of the displacement X{t). We need this in order to exclude 
the possibility that the a priori martingale part of the displacement and the martingale 
approximation of the compensator in the limit just cancel out. (As it is well known, 
this happens for example in tagged particle diffusion in Id simple symmetric exclusion 
process with nearest neighbour jumps, [2].) In paragraph 5.1.2, we prove the i/_i-bound 
(123) for our particular case. Finally, in subsection 5.2, we check conditions (122) of 
Theorem SVY for our particular model. 

5.1 Diffusive bounds 
5.1.1 Lower bound 

For s, t e M with s <t, let 

M{s, t) := X{t) - X{s) - (p{v{u)) du = B{t) - B{s). (124) 

Lemma 1. (1) Fix s e M. The process [s, oo) 3 t ^ M{s,t) is a forward martingale 
with respect to the forward filtration {J^(-oo,t] '■ t > s} of the process t ^ vi't)- 
(2) Fix t G M. The process (— oo,t] 9 s h-> M(s, t) is a backward martingale with respect 
to the backward filtration {T[s^oo) '■ s <t} of the process t h- > rjit). 

Proof. There is nothing to prove about the first statement: the integral on the right-hand 
side of (124) was chosen exactly so that it compensates the conditional expectation of 
the infinitesimal increments of X[t). 

We turn to the second statement, which does need a proof. This consists of the 
following ingredients: 

(1) The displacements are reverted on the flipped backward trajectories 1 1— > 'fj{t) deflned 
in (105): 

X{t)- X{s)^ -X{t)+X{s). (125) 

(2) The forward process t r]{t) and flipped backward process t h- >■ rj{t) are identical in 
law (Yaglom reversibility). 

(3) The function uj i— > (p{u}) is odd with respect to the flip-map uj i— > —uj. 
Putting these facts together (in this order), we obtain 

lim /i-iE(X(s -h)- X{s) I J^[,,oo)) = lim ^"^E( - X(-s + h) + X{-s) \ ^(-oo,~s]) 

= -^ivi-s)) = ifiivis)). (126) 

□ 
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From Lemma 1, it follows directly that for any s < t, the random variables M{s,t) 
and (p{'r]{u)) du are uncorrelated, and therefore 

E{{X{t) - X{s)r) = E{{M{s, t)r) ^iviu)) duY) (127) 

^{t-s) + E{{j\{r,{u))duf). 

Hence, the lower bound in (40). 
5.1.2 Upper bound: if_i-bound 

We recall a general result proved in [19]. See also the surveys [16], [11] and further 
references cited therein. 

Let i I— > C(0 be the reversible Markov process on the same state space (f2, tt) as the 
original rjlt) which has the infinitesimal generator —S. 

Lemma (SVY). Let ip e C^iVL^Ti) with /<^d7r = 0. Then 

limsupriE(( / Lp{ri{s)) dsf) < lim riE(( / ip{^{s))dsf) = 2\\S-^/^ipf . (128) 

t-*oo Jo Jo 

In our case, 

S^-Ia, (129) 

and the reversible process t ^ ^{t) is the so-called diffusion in random scenery process. 
That means: 

m ■= Tz,^ (130) 

where t t— is a Brownian motion in M'^ of covariance independent of the field lo. 
The function : — > M is <^(a;) = <^(0). Thus, the upper bound in (128) will be 

limt-^E(( / ip{i{s))dsf\ = \Ymt-^E{{ I u{Z,)dsf\=l |pP^V(p)dp. (131) 
Jo Jo JRd 

Here, the last step is straightforward computation with expectation taken over the Brow- 
nian motion Z(t) and over the random scenery u. The integral on the right-hand side is 
the same as in (41), and thus, (131) yields the upper bound in (40). □ 

5.2 Graded sector condition 

As a first remark, note that condition (122) is equivalent to 

II S-'/''A+S''/^ U II < Cn\ (132) 
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where the operator S ^^^A+S is meant as first defined on a dense subspace of 

Tin and extended by continuity. In our case, the dense subspace will be C„ specified in 
(89) and 



-1/2 



(133) 



The operators | A | map the subspaces C„ to themselves and are bounded, see (79), 
(80). In order to bound the norm of the operator | A | a\ : Tin Ti-n+i, let u & Cn, 
then 

i 1 

I A p^/^ aiU(pi, . . . ,Pn+l) = . — , , ^ r ^^^^^^^(Pl, ■ ■ ■,Pm-l,Pm+l, ■■■,Pn)- 

(134) 

(135) 



Hence 
(n+1) 



1-1/2 *^ 
I CiiU 

1 



n+l 



n+l 



n 



=1 I 



dpi...d]?. 



< in 



= (n + l) 



y^Pm;^(Plv,Pm-l,Pm+lv,Pn) 
m=l 

+ I ^l^^n+Ml^l^(^^l'->Pn)|^ n ^^^^dpi...dp„ 

U'i JR'i lE^tlPml m=l 

/ f l-f M2 TT ^(P"*) I f ^"+1'' ^(Pn+l) 
/ •••/ \u{pir-,Pn) I 11 ^ -2 

JtR" JR<i ^=1 |Pm| 



n+l 



+1 „ |2 



dpn+1 dpi...dp„. 



In the second line, Schwarz's inequality and the symmetry of the function u{pi, . . . ,Pn) 
is used. 

The innermost integral of the last expression in (135) is bounded above by 

V{p + q) 



C := sup 



Q\ 



dq < oo. 



Thus, for u & Cn 



-1/2 



a,u 



< C^(n+1) \\u 



Hence, by continuous extension. 



and (132) with 7 = 1/2 follows. 



-1/2 



«* l-Hr, 



< CVn + 1, 



(136) 

(137) 

(138) 
□ 
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